Dust in the inner regions of debris disks around A stars 
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ABSTRACT 

We present infrared interferometric observations of the inner regions of two A-star debris 
disks, j3 Leo and £ Lep, using the FLUOR instrument at the CHARA interferometer on both 
short (30 m) and long (>200 m) baselines. For the target stars, the short baseline visibilities 
are lower than expected for the stellar photosphere alone, while those of a check star, S Leo, 
are not. We interpret this visibility offset of a few percent as a near-infrared excess arising from 
dust grains which, due to the instrumental field of view, must be located within several AU of 
the central star. For (3 Leo, the near-infrared excess producing grains are spatially distinct from 
the dust which produces the previously known mid-infrared excess. For £ Lep, the near-infrared 
excess may be spatially associated with the mid-infrared excess producing material. We present 
simple geometric models which are consistent with the near and mid-infrared excess and show 
that for both objects, the near-infrared producing material is most consistent with a thin ring 
of dust near the sublimation radius with typical grain sizes smaller than the nominal radiation 
pressure blowout radius. Finally, we discuss possible origins of the near-infrared emitting dust in 
the context of debris disk evolution models. 

Subject headings: circumstellar matter — stars: individual(beta Leo, zeta Lep) 



1. Introduction 

The list of main sequence stars known to have 
circumstellar material in the form of debris disks 
has been greatly expanded over the last few years 
by surveys at longer wavelengths and most re- 
cently fr om Spitzer observa tions (see e.g., the re- 
view by iMever et al.l 12007). Given the size and 
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distribution of dust in these disks, the grains are 
expected to have short lifetimes. Therefore, it is 
generally believed that the dust is not remnant 
from the star formation process, but is generated 
through collisions of larger bodies. The majority 
of known debris disks have cold (<100 K) mate- 
rial located tens of AU from the central star in 
an analog of our own Kuiper belt. In some cases, 
this material extends to 1000 AU. A small frac- 



tion ( Rieke et a l. 2005; Beic hman et al . 200G) have 
warmer dust located within 10 AU of the central 
star. 

The distribution of material in a debris disk is a 
balance of collisions, radiation pressure, Poynting- 
Robertson (PR) drag and the dynamical influence 
of any large bodies in the system. In order to con- 
strain models of these systems, the dust spatial 
extent and grain size distribution must be mea- 
sured. Observations of optical and near-infrared 
scattered light have provided the most detailed 
overall picture of the dust distribution. However, 



1 



these scattered light observations do not have suffi- 
cient resolution to characterize the material closest 
to the star, and this is where infrared interferom- 
etry can provide a unique constraint. 

Although many of these sources do not show a 
clear near-infrared excess in their spectral energy 
distribution (SED), limits set by spatially unre- 
solved broadband photometry are generally not 
better than a few to several percent. A small, 
warm dust component could be present if dust 
generated by collisions migrated close to the star 
or was produced by bodies in close orbits. If lo- 
cated within a few AU of the central star, this dust 
would be at temperatures which would produce 
near-infrared emission and small grains would pro- 
duce scattered light. Detection of (or stringent 
limits on) warm dust will characterize the inner 
portions of these debris disks. The spatial resolu- 
tion of infrared interferometry can be exploited to 
probe for warm dust in these systems. On long 
baselines (> 100 meters) the central star is re- 
solved and the visibility is primarily a measure 
of the stellar photospheric size. On shorter base- 
lines (< 50 meters) the photosphere is mostly un- 
resolved and if the measured visibilities have high 
accuracy, other emission components can be de- 
tected by looking for deviations from the visibility 
expected for the stellar photosphere. Any resolved 
or incoherent emission will decrease the measured 
visibility from the stellar value. 

Teams using th e Palomar Testbed Interfer- 



ometer ( Ciardi et al.l 2001) and the Center for 



High Angular Resolution Array ( Absil et al. 2006t 



di Folco ct al. 2007; Absil et al 



2008!) have de- 
tected near-infrared extended emission around 
known debris disk systems, including Vega, the 
prototype debris disk. While a near-infrared ex- 
cess was not known through broadband spectral 
modeling, the interferometrically detected near- 
infrared excess was consistent with the photomet- 
ric uncertainties. Observations of other debris 
disk sources revealed a small near-infrared excess 



flux around r Ce ti and C Aql (|di Folco et al.ll2007 ; 
Absil et al.l l2008h . In all these systems, the near- 
infrared excess is consistent with emission from an 
inner, hot dust component, although for £ Aql, a 
binary companion is also a likely origin. 

In this paper we present infrared interferometry 
observations of two known debris disk systems, the 
A- type stars, (3 Leo and £ Lep. The interferom- 



etry observations, including determination of the 
stellar angular diameter, and mid-infrared imag- 
ing for (3 Leo are presented in <j2] Possible origins 
for the observed visibility deficit are discussed in 
SJ21 In 21 we discuss the distribution of the ex- 
cess producing grains and in §5 the origin of these 
grains. Our conclusions are given in §6. 

2. Observations and data analysis 
2.1. Targets 

The targets were chosen from the sample of 
known debris disk systems with the V and K 
brightness as the primary selection criteria. Table 
[TJ lists the target and check star stellar properties. 

(3 Leo was identified as hav ing an infrared ex- 
cess from IRAS observations ( Aumann fe Probstl 
199lh . Mid-infrared imaging ha s not resolved 



the disk ( Javawardhana et al ]l200ll . gUt although 
differences between t he IRAS and ISO fluxes led 
Laureiis et al. (2002) to suggest that the disk 



emission may be so mewhat extended in the ISO 
beam (52" aperture). IChen et all (|2006l ) obtained 
Spitzer IRS spectra of (3 Leo and found a fea- 
tureless continuum spectra consistent with dust 
at ~120 K located 19 AU from the central star. 
C Lep was also i dentifi ed as a debris disk by 



Aumann fe Probstl (|199lh and has an unusually 



high d ust temperature (>30 K) ( Aumann fe Probstl 



19911: IChen fc Jural l200lh . Recent 
Moerchen et alj (|2007l ) resolved the 



work by 

excess at 



18 /im and their model comprises two dust bands 
extending from 2 to 8 AU. As with j3 Leo, the 
Spitzer IRS spec trum for C, Lep is featureless 
(jChen et al.ll200fih . 



2.2. CHARA Observations 

Observations were conducted with the FLUOR 
fiber-optics beam combiner at the Center for 
High Angular Resolution Array (CHARA) op- 
erated by Georgia State University. CHARA is a 
long-baseline, six-element interferometer with di- 
rect detection instruments that work at optical to 



near- infrared wavelengths (jten Brummelaar et al 
2005). These FLUOR observations were taken 
in the K' band and have an effective central 
wavelength of 2.14 microns. The FLUOR beam 
combiner produces high precision visibilities by 
interfering the inputs from two telescopes af- 
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Parameter 


(3 Leo 


S Leo 


C Lep 


HD number 


102647 


97603 


38678 


Spectral type 1 


A3Va 


A5IVn 


A2Vann 


Distance 2 (pc) 


11.1 ±0.11 


17.7± 0.26 


21.5 ±0.32 


Radius 3 (R ) 


1.58 ± 0.018 


2.17 ± 0.073 


1.60 ± 0.11 


T c ff 


9020 4 


8296 1 


9910 1 


Luminosity 5 (L©) 


11.5 ± 1.1 


15.5 ± 1.8 


17.0 ± 2.3 


v sin i (km/sec) 


no 4 


173 6 


245 4 



Table 1 : Stellar properties of the sources 



X NASA Stars and Exoplanet database: http://nsted.ipac.caltech.edu 

2 Dista nces taken from Hipparcos llPerrvman &; ESA 1997); we note a more recent reduction of Hipparocs data l lvan Leeuwenl 
2008) has yielded new distances which are within 1 a of those listed here. We use the older values for consistency with previous 
work. 

3 this work 

1Chen et all ll2006h 

5 calculated from the radius and effective temperature 
qRieke et all (120051 ') 



ter spatial filtering through single-mode fibers 
( Coude du Foresto et al . 2003). In this design, 
wavefront aberrations are converted to photomet- 
ric fluctuations which are corrected by simultane- 
ous measurement of the fringe and photometric 
signals from each telescope. 

P Leo and a check star S Leo were observed on 3 
nights in 2006 May and £ Lep on 2 nights in 2006 
October and November. A check star for £ Lep 
was observed but due to its lower K band flux, 
these data were not useful and are not included 
here. Observations of the targets and check star 
were interleaved with calibration observations to 
determine the instrument response function, also 
called the system visibility. The check star is an 
additional target with roughly the same properties 
as the main target, but no known excess emission 
at any wavelength. Observations of the check star 
are processed in the same way and with the same 
calibrators as the main target and serve as a mea- 
sure of systematic effects in the data. The calibra- 
tors used, along with their adopted diameters are 
given in Table [2] 

The FLUOR data consist of temporally modu- 
lated fringes over an optical path difference (OPD) 
of 170 microns, centered around the zero OPD. 
The coherence length (fringe packet size) in the 
K' band is of order 11 fringes, or approximately 25 
microns. In addition to the fringe signal, FLUOR 



records simultaneous photometric channels, in or- 
der to allow the correction of scintillation noise 
and coupling variations in the input single mode 
fibers. The photometric correction and flux nor- 
malization were done using the numerica l meth - 



ods described in ICoude Du Foresto et al.1 (|1997D . 



Once the fringe signal was recovered, we estimated 
the squared visibilities of individual frames as the 
integrated power in the frequency domain. 

To estimate the fringe power, we used a 
time/frequency transform, a Morlet wavelets 
transform, i nstead of the classical F ourier ap- 



proach (jCoude Du Foresto et al.lll997h . The clas- 
sical Fourier method extrapolates the power under 
the fringe peak using d ata collected at frequ encies 
outside the fringe peak ( Merand et al.|[2006h . This 
approach works well if the readout noise is white. 
The wavelets approach isolates the fringe signal 
in the OP D and in the frequenc y domains (as de- 
scribed in lKervella et al.ll2004al ). allowing a mea- 
sure of the off-fringe power at all frequencies and 
therefore a direct measurement of the background 
noise for each scan. The isolation of the fringe 
signal in OPD is possible because the modulation 
length used (170 microns) is much larger than the 
coherence length (approximately 25 microns) and 
the background noise is measured using the por- 
tion of the scan situated more than 50 microns on 
each side of the fringe packet (i.e. four times the 
coherence length). 



3 



Calibrator Diameter (mas) Target 



Diameter reference 



70 Leo 
C Vir 

IRC 10069 
77 Lep 
HR 1965 
HR 1232 



0.770 ± 0.015 
0.760 ± 0.015 
1.342 ± 0.07 
0.940 ± 0.020 
1.272 ± 0.017 
0.920 ± 0.020 



/3 Leo, 
(3 Leo, 
C Lep 
C Lep 
C Lep 
C Lep 



5 Leo 
5 Leo 



SB relation, Kervella ct al. 

SB relation. Kervella et al. 

Merand et al. (2005) 

SB relation. Kervella et al. 

Merand et al. (2005) 

SB relation, Kervella et al. 



(2004b) 
(2004b) 

(2004b) 

(2004b) 



Table 2: The calibrators used for the CHARA observations. The cali brator sizes are derive d using optical 
and infrared photom etry and the surface brightness (SB) relation from lKervella et ah ( 2004b ') or taken from 
Merand et ail (|2005l ). 



The background noise arises from 3 compo- 
nents: the photometric variation residuals (after 
photometric correction), the photon shot noise 
and the detector readout noise. The first com- 
ponent is only present at very low frequencies, 
since fringes are acquired at a frequency (100Hz) 
higher than the scintillation and coupling varia- 
tions (typically 25Hz at CHARA) and because the 
photometric correction is very efficient. The sec- 
ond component (photon shot noise) is white noise. 
The third component, readout noise, is less pre- 
dictable and can have transients or peaks at dis- 
crete frequencies (electronic noise) . As the wavelet 
approach directly measures the background com- 
ponent from the data, there are fewer residuals 
than in the Fourier method where the noise esti- 
mate is approximate. For the FLUOR data, the 
wavelet method improved the consistency of the 
results, although the basic results are the same 
between the two methods. 

Finally, the final squared visibility estimate and 
the one sigma uncertainty for a given batch of 
frames are obtained by the average and stan- 
dard devi ation of the boo t strapp ed average, as de- 
scribed in Kervella et alj ( 2004a ). The calibrated 
target data obtained using this reduction method 
are given in Table [3l 

2.3. Stellar size and visibility deficit 

If the measured visibilities were due entirely to 
a resolved stellar disk, both the short and long 
baseline data would be well-fit with a single uni- 
form disk. However, as shown in Figure[TJ the visi- 
bility measured on the short baseline for (3 Leo and 
C, Lep is lower than expected from the stellar size 
fit on the long baseline. Fitting a single stellar size 
to both baselines yields a very poor fit as measured 



by \r m comparison to the single-baseline only fits 
for the target stars, while the single-component fit 
to both baselines for the check star, 5 Leo, is good 
(Tabled]). Any additional flux component within 
the field of view will decrease the measured vis- 
ibility and will therefore make the model more 
consistent with the short-baseline data. A par- 
tially resolved emission component will increase 
the discrepancy between the long and short base- 
line visibilities as it would be more resolved, and 
therefore have lower visibility, on the long base- 
lines. An over resolved, i.e. incoherent, source of 
emission will produce the same fractional decrease 
in visibility for all baselines. For the simple case of 
a star and an incoherent component, the measured 
visibility, V mcas , is 



V, 



2 

mcas 



Kstar * /si 
/star / 



incoh 



(1) 



where T4tar is the visibility of the stellar photo- 
sphere and /star arid /i llco h are the fractional stel- 
lar and incoherent component fluxes. The visibil- 
ity used here is a normalized quantity such that 
an unresolved source has V — 1 while an incoher- 
ent (i.e. completely resolved) source has V — 0. 
We fit a single uniform diameter plus an incoher- 
ent emission contribution to both baselines, which 
gives a lower \r f° r P Leo and Q Lep than the 
uniform disk by itself. The visibility deficit for 
( Lep is a tentative detection as the V 2 predicted 
from the stellar size, 0.996±0.001, is only 2.3 a 
from the average measured visibility, 0.966 ± 0.013 
and the stellar size uncertainty is much larger as 
the star is smaller and fainter than the other tar- 
gets. We note that the stellar size for £ Lep from 
data reduced using the classical Fourier approach 
is the same as for the wavelets approach, despite 
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Object 


MJD 


Baselinc(m) 


Pos Angle (deg) 




er 


f3 Leo 


53856.226 


32.531 


-12.980 


0.9487 


0.0219 




53856.270 


33.234 


-21.034 


0.9001 


0.0269 




53856.309 


33.801 


-26.735 


0.9285 


0.0204 




53864.185 


313.083 


74.485 


0.0679 


0.0079 




53865.185 


312.858 


74.288 


0.0503 


0.0040 




53865.236 


293.126 


68.316 


0.0897 


0.0032 


(5 Leo 


53856.248 


33.823 


-20.956 


0.9726 


0.0393 




53856.290 


34.069 


-27.244 


1.0025 


0.0173 




53856.328 


33.940 


-31.660 


1.0356 


0.0482 




53864.233 


286.055 


62.950 


0.2206 


0.0253 




53865.215 


295.815 


65.870 


0.1790 


0.0088 


C Lep 


54040.479 


218.336 


-57.363 


0.7543 


0.0713 




54040.487 


223.077 


-57.415 


0.9403 


0.1002 




54040.506 


232.879 


-57.237 


0.5783 


0.0772 




54045.475 


24.739 


-22.258 


0.9742 


0.0209 




54045.498 


26.062 


-26.679 


0.9524 


0.0235 




54045.518 


27.252 


-29.651 


0.9696 


0.0272 



Table 3: The calibrated visibility observations from CHARA. 



the scatter in the long baseline data. Further ob- 
servations are needed for confirmation of the vis- 
ibility deficit of C Lep. The best- fit incoherent 
component corresponds to an excess flux of 2.7 ± 
1.4 Jy for [3 Leo and 0.47 ± 0.41 Jy for C Lep. For 
S Leo the uniform disk fit is adequate, suggest- 
ing no visibility deficit on the check star and no 
substantial systematics in the observing or data 
reduction process. 

From the uniform disk fit, we can calcu- 
late a limb-darkened angular and physical di- 
ameter for these star s usin g the formula from 
(119741 ) 



Hanburv Brown et al 



"ld 



1 - Ma/3 
1 - 7 M a/15 



1/2 



(2) 



where the coefficient /i depends o n the effec- 



tive t emperature and is taken from iClaret et al 
(jl995l ). The difference between the limb darkened 
and uniform disk diameters is less than 2% for 
our stars, with ratio values ranging from 1.011 
to 1.014. The uniform disk and limb-darkened 
diameters and the derived stellar radii are given 
in Table [51 where the uniform diameter is taken 
from the stellar + incoherent component model. 
These limb-darkened diameters agree with values 
calculated from the su rface-brightness relation of 
Kervella et al] (|2004af) of 1.35, 1.17 and 0.73 mas 



for (3 Leo, 6 Leo, and ( Lep respectively. 

We note that these limb-darkening parameters 
are appropriate for slowly rotating stars, which is 
vio lated by the values for y sin i given in Table 
[TJ lAufdenberg et all ( 2006f ) find limb-darkening 
corrections 2.5 times higher for Vega, an AO star 
rotating at 275 km/sec. As our data are insuf- 
ficient to separately derive the limb-darkening or 
rotational velocity, we use the low rotation rate co- 
efficients to allow for comparison to other works, 
but note that even at 2.5 times higher, the limb 
darkening corrections would be 4%, still much too 
small to explain the difference between the short 
and long baseline sizes given in Table [4] 

The diameter of (3 Leo has been previ ously mea- 



sured with intcrferometry observations. lHanburv Brown et al 
([19741 ) obtained a limb-darkened diameter of 
1.33 ± 0.1 mas at a wavelength of 4430 A 
with the Narrabri in tensity interferometer, while 
di Folco et all (|2004 measured 1.449 ± 0.027 mas 



at 2.17 [im with the VLTI, which is inconsistent 
with our diam e ter at the 3.7<r level. However, the 



di Folco et al.l (|2004h fit did not in clude an inco 



herent component. If we include the ldi Folco et al 



(|2004|) data in our two component fit, both the 
stellar diameter and incoherent flux level change 
by less than 0.2 a, thus the VLTI and CHARA 
data are consistent. 
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Object 


Uniform disk 




Uniform disk 




Uniform disk 


Uniform disk - 


- incoherent flux 






all data 




long baselines 


short baselines 


all data 






Diam.(mas) 


Xr 


Diam.(mas) 




Diam.(mas) 


x? 


Diam.(mas) 


Inc. flux 


Xr 


(3 Leo 


1.332 ± 0.014 


3.8 


1.332 ± 0.009 


2.1 


2.289 ± 0.31 


0.9 


1.323 ± 0.013 


0.024±0.013 


1.6 


6 Leo 


1.148±0.025 


0.8 


1.149 ±0.012 


0.5 


0.0 ± 1.17 


0.5 


1.149±0.022 


0.0±0.006 


1.0 


C Lep 


0.70±0.15 


2.5 


0.69 ±0.09 


3.7 


2.0 ± 0.65 


0.2 


0.66 ±0.14 


0.015±0.013 


2.0 



Tabic 4: Uniform diameter and incoherent flux fit to data 
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Fig. f. — The measured visibilities and errors (points) for (3 Leo (top, left), S Leo (top, right) and 
£ Lep(bottom, left). For each object, the visibility curves for a uniform disk fit only to the long base- 
line data (solid line with dotted line errors) and for a uniform disk plus an incoherent flux component fit to 
all the data (dashed line with dot-dash line errors) are shown. 



Object 


Uniform disk diam. 


LD coeff 


Limb darkened diam. 


Stellar radius 




(mas) 




(mas) 


(R©) 


f3 Leo 


1.323 ± 0.013 


1.012 


1.339 ± 0.013 


1.54 ± 0.021 


6 Leo 


1.149 ± 0.022 


1.014 


1.165 ± 0.022 


2.14 ± 0.040 


C Lep 


0.66 ± 0.14 


1.011 


0.67 ± 0.14 


1.50 ± 0.31 



Table 5: Measured uniform disk and limb-darkened diameters 
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2.4. Mid-infrared imaging 

Mid-infrared imaging observations of [3 Leo 
were made on 2006 March 8 (UT) using the 
Mid-Infrared Echelle Sp e ctrom eter (MICHELLE; 



1993) on the Gemini 



Glasse fc Atad-Ettedgui 
North 8-meter telescope. MICHELLE utilizes a 
320 x 240 pixel Si:As blocked impurity band de- 
tector, with a spatial scale of O'.'l pixel -1 . Imag- 
ing was obtained in the Qa filter (A c = 18.1/xm, 
AA = 1.9/xm) with a standard off-chip 15" ABBA 
chop-nod sequence and a chop position angle of 
30 deg E of N. Two image sequences of (3 Leo were 
taken with 30 ms frametimes and a total on-source 
integration time of 325 s per image. Prior to and 
following the (3 Leo observations, HD109511 (K0, 
-Pi8/im ~ 1-4 Jy) was observed with the same ob- 
serving sequence to serve as a point spread func- 
tion and flux density calibrator. The data were 
reduced with custom-written IDL routines for the 
MICHELLE data format. 

f3 Leo appears unresolved in comparison to the 
calibrato r. At 18.5 p, t he excess for (3 Leo is 
~ 0.3 Jy (|Chen et alJl2006h and at 19 AU from the 
star, the radius inferred bv lChen et al. I (|2006h for 
the mid-infrared emitting material, we measured 
an rms dispersion of the background in the j3 Leo 
images of ~ 0.55 mJy pixel -1 . In §4.41 we will 
use this limit to constrain the radial extent of the 
mid-infrared emitting material. 

3. Origin of the visibility deficit 

In this section we discuss the possible origins of 
the visibility deficit. 

3.1. Companion 

A companion anywhere within the 0'.'8 (FWHM) 
field-of-view (FOV) will lower the measured vis- 
ibility. A companion within the fringe envelope 
(roughly 25 milliarcsec for these observations) will 
produce a visibility modulation which is a func- 
tion of the binary flux ratio and separation and 
the projected baseline length and position angle. 
A companion outside this separation range but 
within the field-of-view will contribute incoher- 
ent flux and the visibility decrease will be the 
same fraction on all baselines. The flux ratio of 
a companion which would produce the measured 
visibility is the incoherent fraction listed in Table 
IH which corresponds to AK = 4.0 ± 0.9 for (3 Leo 



and AK = 4.5 ± 1.4 for ( Lep. These flux differ- 
ences would be produced by a main sequence star 
of spectral type M0 for (3 Leo and M2 for ( Lep. 

Neither star has a known companion within a 
few arcsec of the primary star. The Washington 
Double Star (WDS) catalog lists 3 companions 
for (3 Leo, located from 40" to 240" from the 
primary (far outside the F OV) with V magni- 
tude differences of 6.3 to 13 ( Worlev fe Douglass] 
19971 ). None of these stars could affect the in- 



terferometry observations due to the large an- 
gular separation. Q Lep has no listed compan- 
ions in the WDS. Bot h objects have been imaged 
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in the mid-infrared (jJavawardhana et al 
Moerchen et al. 1 120071 . g23) with no companion 



detected. In our MICHELLE/GEMINI data, the 
Q-band magnitude difference for a point source 
which can be ruled out is 2.5 mag within 0"5 and 
4 mag from 0"5 to 0"8. These data are sufficient 
to detect a possible companion between 0"5 to 0"8 
around {3 Leo for the derived companion spectral 
type of M0. 

The strongest constraints on close (<1") com- 
panions come from the Hipparcos measurements. 
(3 Leo was observed 64 times over 3.0 yrs with fi- 
nal positional uncertainties of 0.99 mas (RA) and 
0.52 mas (dec), and ( Lep was observed 117 times 
over 3.1 yrs with un certainties of 0.51 mas (RA) 
and 0.41 mas (dec) (jPerrvman fc ESAlll997h . As 
neither source was detected to have any astromet- 
ric motion by Hipparcos, these uncertainties can 
be used to place limits on any stellar companions. 
Using the secondary stellar types inferred from the 
flux ratios, the companion stellar masses would be 
approximately 0.5 M Q for (3 Leo and 0.4 M Q for 
C Lep. As the astrometric signature increases with 
orbital distance, the astrometric uncertainty from 
the Hipparcos data sets a lower limit to the ex- 
cluded periods, while the sampling duration sets 
the limit for longer-period companions. To es- 
timate the shortest period companion which the 
Hipparcos data could detect, we assumed a mass 
for each primary of 2.0 M Q and quadratically com- 
bined the positional uncertainties to obtain astro- 
metric uncertainties of 1.12 mas for (3 Leo and 0.65 
mas for ( Lep. Setting a threshold of 5er to ac- 
count for the uneven time sampling, the minimum 
detectable separations are 0.25 AU ((3 Leo) and 
0.34 AU (C Lep), which correspond to periods of 
32 days and 51 days respectively. 
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The detection of longer-period companions is 
limited by the overall time span of the Hipparcos 
data. The orbital period and astrometric signa- 
ture of a companion located at the edge of the 
FOV would be 5.5 yrs and 250 mas for (3 Leo and 
14.5 yrs and 200 mas for £ Lep. For j3 Leo, the 
Hipparcos data samples half a period and would be 
sufficient to detect such a companion. For £ Lep, 
the Hipparcos data would sample 20% of the or- 
bital period. For a circular orbit, the deviation of 
this arc from a best-fit straight line would be 12 
mas, detectable with the 0.65 mas uncertainty, but 
detecting some phases of an elliptical orbit would 
be more difficult. A very long period companion 
with the relevant magnitude difference could have 
escaped detection if the orbit is inclined on the 
sky such that companion is currently too close to 
the primary (within 0'.'5) for detection by imag- 
ing. One probe of such a very long period or- 
bit is the proper motio n as a function of time. 
Gontcharov et al.l (|200l[ ) combined proper motion 
data from ground-based catalogs starting in the 
1930's with the Hipparcos data. For (3 Leo and 
£ Lep, the combined proper motions were within 
the uncertainties of the Hipparcos proper motions, 
and both stars were classified as having no com- 
panions within 10". 

Any companion closer than the short-period 
limit derived above would produce a substantial 
radial velocity signature. Using the inclination an- 
gles derived in 93.21 a companion at the short pe- 
riod limits above would produce a radial velocity 
of 8 km/sec for (3 L eo and 12 km/sec for £ Lep. 
Galland et al.l (|2005l ) made measurements of Leo 
with an uncertainty of 137 m/s, more than suffi- 
cient to detection such a large signature, however, 
the time sampling covered only a few hours and 
is not sufficient to rule out compan ion periods of 



tens of days. Observations of £ Lep (jGrenier et al 
19991 , e.g.) have also been made with sufficient 



precision, but not sufficient time sampling to find 
a companion with a period of a many days. 

In summary, neither target star has a known 
companion within the CHARA FOV and Hippar- 
cos measurements rule out companions with pe- 
riods from tens of days to several years. A very 
close companion (periods less than tens of days, 
separations less than 0.35 AU) can not be ruled 
out in either case, but would produce an easily 
detectable (> 5 km/sec) radial velocity signature. 



Although we can not definitively rule out a com- 
panion as the source of the flux decrement, it is 
unlikely given the above constraints on period and 
magnitude difference. A less massive companion 
would produce a smaller flux decrement, which 
would require another flux component in the sys- 
tem. Given the small phase space remaining for an 
undetected companion and the fact that the two 
mid-infrared excess sources ([3 Leo and £ Lep) have 
a near-infrared visibility decrement, while 6 Leo 
with no mid- infrared excess does not, we proceed 
with the hypothesis that the flux decrement does 
not arise from a companion. 

3.2. Stellar rotational oblateness 

Our analysis of the visibility deficit on the short 
baseline relies on knowledge of the stellar size from 
the longer baselines. If the star is oblate due to 
rotation, the predicted size on the short baseline 
may be incorrect as the short and long baselines 
are nearly orthogonal (Table [3]). We can calcu- 
late the maximum possible effect by assuming the 
short stellar axis is aligned with the longer base- 
line, which would place the longer axis along the 
short baseline, producing lower visibilities. We 
calculate the ratio of stellar radii, Xr from 



X R = 



R 



pol 



FL 



V cq R C q 

2GM 



(3) 



where R po i and R e „ are the polar and equato- 



rial radii, 



is the equatorial velocity, G is 



the gravi tational constant and M is th e stel- 
lar mass ( Domiciano de Souza et all l2002h . For 



the most conservative calculation, we take v eq to 
be the max i mum equatorial velocity inferred by 
Rover et al.l ( 2007 ) of a survey of A stars, which 
are grouped by sub-class. These velocities are 300 
km sec -1 for (3 Leo and £ Lep and 280 km sec -1 
for 8 Leo. The resulting oblateness is corrected 
for viewing angle by deriving i from the measured 
v sin i and the assumed v cq and approximating the 
stellar shape as an ellipsoid (Table [6]). The ob- 
served stellar radii ratio X \, s is then given by 



Xr 



(1-(1-X R ) cos 2 ;) 1 ^ 



(4) 



Starting with the derived stellar size on the long 
baseline (0\ ong , see Table HJ, we calculated the 
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V 2 that would be measured on the short baseline 
(V 2 (#ion g ))- We then applied the observed oblate- 
ness factor, A bs to find the maximum possible 
angular diameter, 9i on g/ X ^ s , and recalculated the 
V 2 for the short baseline (V 2 (0i ong /^obs))- Be- 
cause these angular sizes are at best marginally 
resolved on the short baseline, the change in visi- 
bility is less than 1% in all cases, even if the appar- 
ent angular size changes by 20%, as predicted for 
£ Lep. For comparison, we also list the short base- 
line size, # s hort from Table[4j The measured visibil- 
ity on the short baseline, V 2 ncasurod is significantly 
lower than either V 2 (#i on g) or V 2 (0\ ong /X o b s ) for 
both Leo and ( Lep but not for the check star 
5 Leo and thus stellar oblateness can not account 
for the measured visibility deficit. We note that if 
rotational axis of the star is aligned such that the 
short stellar axis is along the short baseline, then 
the true visibility decrement is actually slightly 
larger than measured. 

As these stars are rotating rapidly, they are 
also subject to gravity darkening, which produces 
a decrease in the effective temperature from the 
pole to the equator. Since the limb-darkening de- 
pends on the effective temperature, this effect is 
also linked to the apparent oblateness. However, 
this effect is very small compared to the oblateness 
derived above. Usin g the effective tempera ture 
difference found by lAufdenberg et aL (|2006h for 
Vega, an AO star, of 2250 K, the limb-darkenening 
correction for the pole is 0.3% larger than correc- 
tion at the equator. This factor goes against the 
rotational oblateness which makes the equatorial 
radius larger and even with the factor of 2.5 for 
a fast rotating star, is insufficient to explain the 
ratios between diameters fit to the long and short 
baselines of 1.72 ± 0.23 for Leo and 2.9 ± 1.0 
for £ Lep. 

3.3. Emission and scattering from dust 

Dust grains within the field of view will pro- 
duce a near-infrared excess through thermal emis- 
sion and scattering. We assume that there is no 
gas in these debris disks and therefore the inner 
radial limit for the debris disk is the dust sub- 
limation radius. For a sublimation temperature 
of 1600 K and assuming large grains in thermal 
equilibrium emitting as blackbodies, the sublima- 
tion radius is 0.12 AU for Leo and 0.14 AU for 
£ Lep. The 2 /an emission will be maximized for 



dust at the sublimation temperature, so a lower 
limit to the excess luminosit y can be estimated 
following iBrvden et alJ (|2006h 



L. 



dust 



I 



hvT} 



(5) 



where h and k are the Planck and Boltzmann con- 
stants. For a temperature of 1600 K, the frac- 
tional dust luminosity is 2.0±1.1 x 10~ 3 for Leo 
and 9.8 ± 8.5 x 1 0~ 4 f or ( Lep. For compar- 
ison, Chen et all ( 20061 ) calculated mid-infrared 
dust luminosities of 2.7 x 10~ 5 and 6.7 x 10~ 5 
for Leo and ( Lep respectively. However, the 
much larger near-infrared luminosity does not re- 
quire substantially more mass than implied by the 
mid-infrared excess since, as the fractional dust 
luminosity represents the fraction of the star as 
seen by the dust, the calculated fractional lumi- 
nosities are highly sensitive to the dust location. 
An estimate of the minimum mass of near- infrared 
emitting grains can be calculated using the frac- 
tional l uminosity and a ssuming efficiently emitting 
grains (|Jura et al.lll995l ). 



167rL dust 2 
M dust > — — par 



(6) 



where p is the density, a is the grain radius and r is 
the distance from the star. A minimum mass can 
be calculated by using the idust values calculated 
above for small dust grains located near the subli- 
mation radius. For a grain radius of a = lfim, r at 
the dust sublimation radius and p ~ 2 gm cm~ 3 , 
the minimum mass of the near-infrared emitting 
material is 5 x 10~ 9 M m for L eo and 2 x 10~ 9 
Mq for £ Lep. IChen et al. (2006) derived a mass 
for the small grains in the mid-infrared producing 
material of 4.2 x 10~ 6 for f3 Leo and 5.6 x 10~ 6 
for C Lep. So although the near-infrared excess 
represents a higher fractional dust luminosity, this 
can be produced by a much smaller mass than the 
mid-infrared ring. 

A ring of hot dust near the sublimation ra- 
dius is not incompatible with the incoherent flux 
model fit in 32.31 as the sublimation radius is large 
enough to be resolved on even the short baseline. 
Given the relative uncertainty in the incoherent 
flux component fit, a component with V 2 < 0.2 
would fit within the uncertainty. For Leo, the 
sublimation radius corresponds to 11 mas and a 
ring of any width at this radius has a V 2 < 0.2 on 
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f3 Leo 


6 Leo 


C Lep 


vsini (km sec 1 ) 


110 


173 


245 


assumed v cq (km sec -1 ) 


300 


280 


300 


i (deg) 


21.5 


38.1 


54.7 


R-eq (R©) 


1.54 


2.14 


1.5 




0.74 


0.70 


0.74 


X bs 


0.95 


0.84 


0.80 


fiong (mas) (Table ij) 


1.332±0.009 


1.149±0.012 


0.69±0.09 


9ion S /X ohs (mas) 


1.401±0.009 


1.368±0.014 


0.826±0.11 


^short 


2.289 ± 0.31 


0.0 ±1.17 


2.0 ± 0.65 


V 2 on short baseline: 








v 2 ((W 


0.976±0.0003 


0.981±0.0004 


0.996±0.001 


V (#long/^obs) 


0.973±0.0003 


0.973±0.0006 


0.994±0.002 


V 2 

measured 


0.938±0.015 


1.001 ± 0.015 


0.966 ± 0.013 



Table 6: The calculated maximum visibility change due to rotational oblateness. The uncertainties in 
#iong/^obs, V 2 (#i ng), and V 2 (#iong/^obs) include the uncertainty in the measured value of #i ng but not 
the uncertainty in X b s , which is unknown. 



all baselines in our observations. For £ Lep, the 
sublimation radius is at 7 mas and any ring wider 
than 1 mas (0.02 AU) produces V 2 < 0.2 on all 
baselines. If the inclination angles are close to the 
values inferred in Table l3~2l these approximations 
are sufficient. Thus, thermal emission from hot 
dust near the sublimation radius could produce 
the measured visibility deficit. 

At larger angular scales than have been in- 
vestigated with the interferometer (> 1" — 10"), 
debris disks are often detected in scattered light 
at optical and near- infr ared wavelengths (e.g., 
AU Mic and Fomalhaut: iKalas. Liu, k, Matthews! 



2004; IKalas. Graham, fc Clampinll2005l ) and scat- 
tering from within the field of view of the interfer- 
ometer (< 0'.'8) could also produce the observed 
visibility deficit. Scattering in the near- infrared 
will dominate emission for grains at several hun- 
dred degrees, depending on the grain size and 
composition. To investigate the scattering from 
warm dust, we used the debris disk simulatoiQ 
described by IWolJ (|2Q06T) which calculates the 
thermal emission and scattering given the dust 
size, composition and distribution. For example, 
small grains uniformly distributed from 1.0 to 4.6 
AU (the (3 Leo FOV radius) will produce the ob- 
served near-infrared excess given a total mass of 
small grains of 1 — 7 x 10~ 5 M®, depending on 



1 http: / / aida28.mpia-hd.mpg.de / ~swolf/ dds / 



the exact size and composition. This is more than 
1000 times larger than the minimum mass needed 
to produce the excess from hot grain emission. 

As there is no known evidence for a companion, 
we contend that thermal emission and scattering 
from dust grains is the most likely origin of the 
near-infrared excess. This is also consistent with 
our finding that the two sources with a measured 
visibility deficit have mid-infrared excess emission 
while the control star, which has no known excess 
does not have a visibility deficit. In the next sec- 
tion, we explore the constraints on these grains 
and discuss possible mechanisms for their origin. 

4. Dust distribution and small grain origin 
4.1. Dust grain sizes 

Both (3 Leo and £ Lep have a substantial mid- 
infrared excess which has a characteristic temper- 
ature much lower than dust which would produce 
a near-infrared thermal excess and is therefore fur- 
ther from the central star. Many auth ors (see e.g. 
Dominik fc Decinl 120031: IWvattl [20051 and refer- 
ences therein) have studied the dynamics of debris 
disks similar to our targets and have found that 
collisions are dominant over PR drag, i.e grains 
collide and become smaller before PR drag sig- 
nificantly decreases the size of their orbits. Ra- 
diation pressure also plays a role as small grains 
are subject to removal from the system. How- 
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ever, clearing of small 



lute. Krivov et al. (200 



rains may not be abso- 
modeled the (3 Pic disk, 
which has a similar spectral type (A6V) and opti- 
cal depth to the systems discussed here and found 
that although grains at and below the canonical 
blowout radius are depleted compared to a purely 
collisional system, a population of small grains 
persists in their model. For our target stars, the 
radiation pressure size limit is ~2 /im, the PR 
drag timescale at 1 AU is 1000 yrs for 10 /im ra- 
dius grains and the collisional timescale for these 
same grains is 80 years (f ollowing the formula of 
Backman fc Parescd 1993 ). 

A second constraint on the dust size is the lack 
of a significant sili cate feature in the IRS spectrum 
for either source ( Chen et al. 200G, see Figure 2), 
although the excess for j3 Leo is not strong enough 
at 10/zm to provide as strong a constraint as for 
£ Lep, which has excess emission at shorter wave- 
lengths. The lack of a silicate emission feature 
requires the grain population to have radii larger 
than a few microns if composed of silicates or to 
be primarily non-silicate. 

4.2. Modeling approach 

As the interferometer data provide only an up- 
per limit to the visibility and therefore a lower 
limit to the size of the near-infrared flux region, 
the strongest spatial constraint from the interfer- 
ometry data is that the dust must be within the 
FOV. However, there is another strong constraint 
from the measured mid-infrared excess of these 
sources. The dust producing the near-infrared ex- 
cess will also produce mid-infrared excess, with the 
exact flux depending of course on the dust temper- 
ature and opacity. 

We now begin to explore various specific mod- 
els for the distribution of dust in these systems, 
and examine whether these models fit within the 
constraints provided by the near and mid-infrared 
data. In all models, we assume optically thin 
emission for the near and mid-infrared emission. 
In this section, we consider the relative contribu- 
tions of scattering and emission to the near and 
mid-infrared excess flux. For the sc attering, we 
have used the debris disk models of IWolJ (|200(j| ) 
to calculate the emission and scattered light flux 
for various grain radius and radial distributions 
and two example grain compositions. We have 
chosen grain compositions which will produce the 



featureless mid-infrared spectrum seen in the IRS 
data and have substantially different emissivity ra- 
tios between the near and mid-infrared. These two 
populations are silicate grains with radii between 
3 and 10 /jm and graphite grains with radii from 
0.1 to 100 /im. In both cases, we use a distri- 
bution of grain radii, n(a) oc a~ 35 appropriate 
for collisionally dominated disks. For these toy 
models, we concentrated on illustrative cases of 
dust radial distributions and did not modify the 
grain radius distribution for the effects of radiation 
pressure. The possible presence of small grains is 
discussed in more detail in §4.41 For each case, 
the disk mass was determined by scaling to match 
the observed near-infrared excess. These masses 
are significantly higher than the minimum mass 
derived in £13.31 as that estimate assumes the flux 
comes only from small, hot grains which produce 
much more near-infrared emission for the same 
mass than a distribution of grain sizes and tem- 
peratures can. In Table [7] we present the results 
for the two grain populations over several radial 
distributions, listing the ratio of emission to scat- 
tering at 2 (Ltm, the excess flux at 10 and 24 /zm 
and the mass in small grains. All models have 
radial density profiles of n(r) oc r~ 15 . 

As expected, emission dominates for grains 
close to the central star (<1 AU), while scat- 
tering dominates for grains farther away. The 
mass in small dust grains necessary to produce the 
near-infrared excess flux is higher for scattering- 
dominated disks than for emission-dominated 
disks. The scattering-dominated cases produce 
too much mid-infrared flux, in some cases by more 
than an orde r of magnitude . The 24 /im excesses 
measured bv lSu et all (|2006h are 0.46±0.01 Jy for 
(3 Leo and 0.53±0.02 Jy for £ Lep and the 10 /im 
excess from the IRS spectra are 0.002± 0.004 Jy for 
Leo and 0.18±0.01 Jy for £ Lep from lChen et alJ 
( 20061 ). The emission-dominated disks also pro- 
duce too much mid-infrared flux, but not by as 
large a factor. As the models which have sub- 
stantial near-infrared emission have mid-infrared 
fluxes close to the observed values, we assume that 
emission is the primary mechanism for producing 
near-infrared flux. In the following sections, we 
will explore other density distributions and mod- 
els to fit both the near-infrared and mid-infrared 
excesses in detail. 
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(3 Leo 


C Lep 


Radial 


Grain 






M sma ll { 


;r 


^cm/ -^sc 






M sma ii { 


;r 


distribution 


type 




(Jy) 


(Jy) 








(Jy) 


(Jy) 


M e 




R S ub-FOV 


silicate 


6.1 


1.7 


1.2 


1 x 10" 


-5 


8.5 


18 


16 


1 x 10" 


-4 




graphite 


8.6 


5.4 


2.8 


2 x 10" 


-6 


10 


6.7 


4.6 


2 x 10" 


-5 


1.0 AU-FOV 


silicate 


0.019 


231 


254 


4 x 10 


-4 


0.02 


261 


400 


4 x 10 


-3 




graphite 


4.9 


2.0 


1.8 


2 x 10 


-5 


4.1 


32 


35 


2 x 10 


-4 


Rsub-1.0 AU 


silicate 


8.6 


10 


3.9 


1 x 10 


-6 


8.8 


10.5 


4.0 


5 x 10" 


-6 




graphite 


9.2 


4.7 


2.0 


1 x 10" 


-6 


12 


4.7 


2.1 


5 x 10" 


-6 



Table 7: The ratio of emission to scattering flux at 2 /zm, the mass in small grains necessary to reproduce 
the observed near-infrared excess and the 10 and 24 /im flux for several disk models. The radius of the FOV 
corresponds to 4.6 AU for (3 Leo and 8.6 AU for £ Lep. 



4.3. Dust grain distributions 

For both stars, we first considered the hypoth- 
esis that the grains producing the near-infrared 
excess were generated by collisions between larger 
bodies in the belt which produces the mid-infrared 
excess. These grains can then be dragged to- 
wards the central star via PR drag and become 
sufficiently heated to emit at near-infrared wave- 
lengths. For a specific theoretical description of a 
disk in which grains created in collisions in the 
planetes imal belt mi grate inward, we used the 
model of lWvattl (|2005i ) , who calculated the steady- 
state optical depth as a function of radius. In this 
model, the disks are collisionally dominated, but 
a small fraction of the dust created by collisions 
in the planetesimal belt migrates inwards due to 
PR drag and is subject to col lisions as it m igrates. 
Assuming a single grain size, IWvattl (|2005l ) found 
the optical depth as a function of radius to be 



Teff(r) 



T e ff(r ) 



l + 4?7o(l- 



(7) 



where r e ff(ro) is the optical depth of the planetes- 
imal belt at ro , the radius of the planetesimal belt 
and 7/0 is a parameter balancing collisions and PR 
drag, which had a value of 2.4 for (3 Leo and 6.7 for 
( Lep. For 770 = 1 the collisional lifetime equals 
the time it takes a grain to migrate to the star. 
We assume optically thin, blackbody grains dis- 
tributed with the optical depth given by eq. [7] and 
starting at the sublimation radius. The value of 
Teff^o) is iterated until the optical depth within 
the CHARA FOV produces the observed near- 
infrared excess. We then calculate how much mid- 



infrared flux would be produced and compare to 
the measured mid-infrared excess. 

For (3 Leo, the mid-infrared excess spectra 
is well fit by a grain temperature T «120 K, 
which implies a d istance from the star of 19 AU 
(|Chen et al.ll2006h . well outside the FOV of our 
observations (4.6 AU). Applying the model in eq. 
[7] with ro = 19 AU and assuming an emissivity 
wavelength dependence of A~ 2 fi j4.4[) produces 
a 2/10 /zm flux ratio of 1.6, while the observed 
ratio, using our d etection and the IRS data of 



Chen et al 



(|2006h is > 140. A shallower grain 
emissivity function with wavelength will produce 
an even larger discrepancy between this model 
and the data. Thus, the grains which produce the 
near-infrared excess can not come from a smoothly 
distributed population generated by collisions in 
the mid- infrared belt. 

Moerchen et al.l ( 2007 ) resolved the 18 /zm 



emission from £ Lep and modeled the distribu- 
tion as arising from two rings with stellar distances 
from 2-4 and 4-8 AU, which is contained within the 
CHARA FOV of 8.6 AU. Using a radius of 4 AU 
in eq. [7] produces a 2/6 /im excess o f 0.8, while 



the observed excess from our data and Chen et al 



(120061) is 18. Thus this model i s not a good fit 



for C Lep either. We note that iMoerchen et al 



(|2007l ) did not resolve the excess emission from 
C Lep at 10 /im and concluded that the dust pro- 
ducing this excess is interior to the resolved 18 /im 
rings. They surmise that the 10 /im emitting dust 
is migrating inward by PR drag from the belts 
resolved at 18 /zm. 
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4.4. Geometric models 

To further examine the constraints which the 
near-infrared and mid-infrared excess place on the 
dust distribution, we use a geometric model of the 
dust distribution. The data are compared to the 
models in a Bayesian approach designed to con- 
strain the range of valid model parameters, rather 
than finding a single best-fit model. The input 
data are: 1) the near-infrared excess within the in- 
terferometer FOV and the visibility limits f or this 
excess, 2) the IRS data from IChen et all (|2006h 
and 3) the spectral energy distribution from 2 to 
100 /mi from the literature, including a 70 /an 
Spitzer-MIPS measurement (K. Stapelfeldt, pri- 
vate communication). We have constructed an 
SED for each star using photometry from SIM- 
BAD in order to determine the excess flux. The 
stellar template was determined by fitting the op- 
tical and near-infrared photometry to a gr id of 
Kurucz-Lejeune models ([Leieune et allll997l ) cov- 
ering the range of effective temperature and sur- 
face gravity values appropriate for the main se- 
quence stellar types of the target stars. Both stars 
are nearby an d have photospher ic colors consistent 



with A v = (|Chen et alj |2006) 



The basic disk model is an optically thin ring of 
dust. To simplify the calculations, we consider the 
dust to be geometrically thin; however, we note 
that to intercept ~1% of the starlight, the dust 
will need to have a finite vertical height. At a 
radial distance of 0.1 AU, this corresponds to a 
height h with h/r — 0.02. Such a vertical height 
is sm aller than a flared primord ial disk at this ra- 
dius (jChiang k Goldreicbl[l997l h/r = 0.09) and 
smaller than the value h/r ~ 0.05 derived for the 
P Pic dust disk at larger (r > 15 AU) radii and 
thus is plausible. 

The excess flux ratio between 2.2 /zm and the 
shortest IRS wavelength can be used to set a limit 
on the wavelength dependence of the grain emis- 
sivity. For f3 Leo, F(2.2 /mi)/F(10 fun) > 140, 
implying an emissivity decreasing at least as fast 
as A -13 if the emission is from a hot black- 
body. For C Lep, F(2.2 /zm)/F(6 /mi) ~ 18, im- 
plying emissivity proportional to at least A~°' 9 
for hot dust. As we wish to examine the range 
of disk physical parameters, including the grain 
size and composition, which can reproduce the 
near and mid-infrared excess emission, we adopt 



a analytic approximation for the grain proper- 
ties to keep the calculation manageable yet self- 
consistent. We have therefore chosen a power-law 
representation of the r adiative efficiency e, fo llow- 
ing work on (3 Pic b y Backman et al. (1992) and 
Backman k Paresce! (|l993l ). For a given gram ra- 
dius a, the absorption and emission efficiency is 
roughly constant, e ~ 1 — albedo, for radiation 
at wavelengths shorter than a critical wavelength 
A and decreases for wave lengths longer than A Q 
( Backman k Parescel [l993f) . The relation between 
the grain radius a and the critical wavelength de- 
pends on the grain composition and shape and 
varies from X /a ~ 2-7T for strongly absorbing 
grains to A Q /a ~ 1 /2tt for w eakly absorbing grains 
( Backman k Paresce 19931) . 

We assume that the disk is optically thin to its 
own radiation, therefore the stellar radiation is the 
only input. Radiation from the early A spectral 
types observed here (T e g > 9000 K) is dominated 
by wavelengths < 1/im. Our input data are at 2.2 
/xm and longer, therefore we can not constrain the 
value of the critical wavelength below 2 /xm. The 
grain emission efficiency decreases at wavelengths 
much larger than the grain radius and given our 
wavelength constraints, we assume that the ex- 
cess is dominated by grains with a critical wave- 
length of > 1/xm and therefore that the absorb- 
ing efficiency is essentially constant. The emission 
efficiency is assumed to follow a power-law such 
that e e = e (\ /X) q . This formulation of the effi- 
ciency does not account for spectral features, but 
as neither object has such features, the approxi- 
mation is appropriate. We investigate two values 
of q: q = 1 which is appropriate for absorbing 
dielectrics and amphorous materials such as sili- 
cate and roughly matches the silicate population 
considered in §4.21 and q = 2 which is appropriate 
for conductive substances such as pure graphite 
or crystallines and represents the graphite popu- 
lation in §4.21 We then derived the temperature 
of the gr ains as a function of radius f rom the star, 
following iBackman k Paresce (|l993l) . 



T(r) = 468 Iv 5 A~ 1//5 r -2/ ' 5 K q = 1 



(8) 



T(r) = 685 li /6 A- 1/3 r- 1/3 K q = 2 (9) 

(10) 

where L» is in L Q , A G is in microns and r, the 
distance to the star is in AU. Assuming a power 
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law for the radial distributio n, the flux in an ring 
is then (jKoerner et al.lll998l ) 



F(r,X) 



MV>^ (ID 



where t To is the optical depth and D is the distance 
to the star from earth. The input parameters to 
our models are the inner disk radius rj n , the disk 
radial extent Ar, the optical depth, t To , the optical 
depth radial exponent, a and the grain character- 
istic wavelength A G . Unless the grain composition 
varies with disk radius, the values of t To and e Q are 
degenerate. We set e Q = 1, thus r To here repre- 
sents the emission optical depth and is only equal 
to the geometric optical depth if the grains have 
an albedo of 0. 

For (3 Leo and ( Lep we were unable to fit both 
the near-infrared and mid-infrared excess with a 
single ring of dust. This is not surprising as the 
2 /im excess is higher than the excess flux at the 
shortest IRS wavelengths, requiring a decrease in 
emissivity at some intermediate radii. The next 
level of sophistication is to add a second ring of 
dust, with each ring following the physical descrip- 
tion given above. 

For each object, a grid of millions of models 
was calculated and compared to the input data. 
The results for each object are a range of param- 
eters consistent with the data, given all possible 
values for the other parameters. We found that 
some parameters were not well constrained by the 
data and others were degenerate, such as the opti- 
cal depth and the disk radial extent. We initially 
assumed that the inner and outer rings had the 
same radial power-law a and the same character- 
istic grain size, A Q . We were unable to fit both 
the near-infrared and mid-infrared data with a two 
ring model if the characteristic grain size was the 
same in each ring. Fitting for two values of a and 
A within a single grid is computationally very ex- 
pensive, so we fit just the IRS data to a single 
ring to constrain the values of a and A for the 
outer ring. For /3 Leo, the value of a ou t C r is not 
well constrained, and we assume a value of -3/2 
as predicted for co llisionally dominated disks (e.g. 
Kenvon fc Bromley 20051 ). 

For each star, we tried to match all the input 
data with the emissivity power-laws of 5 = 1 or 
2. For j3 Leo, the data can not be matched with 
(5=1, not surprising given the flux ratio between 



2 and 10 /im discussed above. For £ Lep, models 
with 5 — 1 can match both the near-infrared and 
mid-infrared data, but these models require that 
the outer ring extend over 20 A U, which is much 



larger than the extent derived bv lMoerchen et al 



(|2007l ) in their imaging. We therefore place a prior 
constraint of Ar outcr < 15 AU. With this con- 
straint, only (5 = 2 models provide adequate fits. 

The range of parameter values which falls 
within a 67% probability range (corresponding 
to ±1(7 for a normal distribution) is given for each 
target in Table [H An example model for each 
object is shown with the SED and IRS data in 
Figure [2l For j3 Leo, the outer ring r- ln values of 
7.5-15 AU a re sm aller than the 19 AU found by 
Chen et al. (|2006h due to the different tempera- 



ture law we used, but the inner and outer ring 
are clearly separated by a gap of several AU. For 
C Lep, the inner and outer ring are at similar radii 
(< few AU) and although the inner ring must have 
significantly higher opacity to produce the near- 
infrared flux, it is possible to fit the data with 
models in wh ich the inner and outer ring overlap. 
Interestingly, iMoerchen et all (|2007l ) also required 
a higher flux ratio in their inner ring (2-4 AU) as 
compared to the outer ring (4-8 AU). 

The optical depth and radial extent of the inner 
dust ring in these models are degenerate param- 
eters as the constraining data are the 2 /^m flux 
and the lack of strong mid-infrared flux. We have 
deliberately limited these models to be optically 
thin, but we note that the near-infrared flux could 
also arise from a ring with a very small radial 
extent which was vertically optically thick. The 
strongest test of the radial extent for the inner ring 
would be to resolve it interferometrically, which re- 
quires observations on shorter baselines than the 
data presented here. For example, a ring around 
(3 Leo at 0.12 AU with a radial extent r/4, would 
have V^ ing > 0.5 at 2 /im on baselines shorter than 
13 meters and Vj? ing < 0.1 on the 30 meter base- 
line we used. High precision measurements would 
still be necessary given the small flux contribution 
from the ring. 

The relationship between the characteristic 
grain size and the physical grain radius depends 
on the grain composition and the distribution of 
sizes. One specific example of a grain material 
which could be approximated by our 5 = 2 emis- 
sivity model is graphite, which is strongly absorb- 
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ing for grain radii > 0.1 /im ( Draine &: Led 1984). 
To compute the physical grain size for graphite 
we take X /a ~ 2ir. The other factor is the dis- 



tribution of grain radii. Following iBackman et al 



(|19921 ). one method of tying A to the physical 
radii is to find the radius which divides the grain 
population into two equal halves of surface area. 
For a distribution of n(a) oc a -3 5 , this radius is 
4 times the minimum radius, a ~ 4a m i n . Putting 
these two factors together, we have a m ; n ~ \ /8ir. 
For our upper limit on the characteristic size in 
the inner ring A G < 2/im, this corresponds to 
Omin < 0.08/im. 

Grains that small are below the nominal radia- 
tion blowout radius for these stars. However, ra- 
diation pressure may not completel y clear all the 
grains from debris disks like these. iKrivov et al 



(|2000h modeled the /3 Pic debris disk, which has a 
similar spectral type (A6 V) and mid-infrared ex- 
cess to our targets. In their model, small grains are 
constantly created through collisions, particularly 
between particles on stable orbits and those being 
blown out of the system. They found that grains 
smaller than a few microns were depleted com- 
pared to the a -3 - 5 distribution of a collisionally 
dominated disk, but that a significant population 
remained. The resulting overall grain population 
could be approximated by a more shallow slope 
in the distribution, for example, fitting the result- 
ing grain radii distribution with a single power-law 
between 0.1 and 100 /im, results in n(a) ~ a -2 8 . 
Thus the inner rings may contain grains smaller 
than 1 /im, the nominal blowout radius. 

Our model for both targets includes a much 
larger characteristic grain size in the outer ring, 
A D ~ 35 — 50/im, which corresponds to a mini- 
mum size of ~ 1 — 2/im for graphite grains. This 
is roughly the radiation blowout size. The conclu- 
sion from our models that the inner ring contains 
substantially smaller grain sizes than the outer 
ring should be confirmed with more detailed grain 
models, but as discussed more below, may suggest 
either different origins for the grains or different 
dynamics. 

The dust sublimation temperature of 1600 K 
used in these models may be plausible for amor- 
phous grains such as those represented by the 
S = 1 model, but is higher than ge nerally used 
for c rystalline grains (e.g. 1250 K; Bauer et al. 
19971 ) as represented by the 6 = 2 model. A 



sublimation temperature of 1250 K results in 
much poorer fits to the data. At these high 
temperatures, micro n-sized grai n lifetimes will be 
short; for example, lLamvl (| 1974! ) found lifetimes 
of less than 10 4 seconds for 1 /im radius grains at 
1500-1600 K. However, once the grains are very 
small, the grain temperature and lifetime may in- 
crease. In a study of grains with radii < 0.01 /im 
heated through interactions wi th a single photon, 
Guhathakurta fc Draind (|l989h found a broad dis- 



tribution of temperatures with excursions as high 
as 2800 K for g raphite grains and 2050 K for 



silicate grains. Guhathakurta fc Draine (1989) 



calculated the sublimation rates for these grains 
including a correction derived from fluctuation 
theory for finite systems which decreases the sub- 
limation rate by ~ 10 4 . The resulting lifetimes for 
grains with radii from several to tens of Angstroms 
are > 10 2 yrs. 

The combination of sublimation, radiation 
pressure and collisions will result in a grain size 
distribution substantially more complicated than 
the simple power-law often used in debris disks 
and assumed here. The result of all these pro- 
cesses may be a population of small hot grains 
which is constantly created through collisions and 
depleted through sublimation and radiation pres- 
sure. Alternatively, the presence of a significant 
number of small grains may imply origin in a tran- 
sient event, as discussed in the next section. Our 
formulation of the grain temperature and emis- 
sivity efficiency does not properly represent very 
small grains and more detailed models than those 
considered here are necessary to determine if the 
temperatures and lifetimes of sub-micron sized 
grains are consistent with a stable grain popula- 
tion which could produce near-infrared flux ob- 
served here. Emission from small, hot grains has 
been invoked to fit th e spectral energy dist ribu- 



tions of debris disks ( Sylvester et al. 1997[ ) and 



the much mo re massive primord ial disks of Herbig 
Ae/Be stars (|Natta et al.lll993l ). 

In the models presented here, the mid-infrared 
excess for (3 Leo is produced by dust grains located 
r w 12 ± 5 AU from the star. At the pixel scale of 
our MICHELLE/Gemini data (O'.'l pixel" 1 , $2~H} 
and the distance of (3 Leo (d = 11.1 pc), the mid- 
infrared dust emission is located of ^8 pixels from 
the central core of the stellar image (FWHM ^5.4 
pixels); thus, the Gemini observations of (3 Leo 
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should easily resolve the mid-infrared emission if 
the surface brightness is high enough. However, 
the Gemini observations show no evidence of de- 
tecting the outer mid-infrared emission. For a 5 
AU ring width, the dust emission is spread out 
among ^400 pixels, for a surface brightness of 
~0.75 mJy pixel -1 . Given the measured disper- 
sion of ~ 0.55 mJy pixel -1 , the signal-to-noise ra- 
tio is then only ~ 1.4 (per pixel). Thus, the non- 
detection of the mid-infrared excess of Leo is 
fully consistent with our dust distribution model, 
even though the spatial resolution was more than 
sufficient to resolve the emission. 

5. Hot dust in debris disks 

Photometric and spectroscopic surveys of de- 
bris disks have revealed other sources with dust 
within several AU of the star and the frequency 
of these systems is a strong function of the age of 
the system. For FGK stars the frequency of debris 
disks emitting at wavelengths shorter than 30/im 
is 9-19% at ages less than 300 M yr, but less than 



2006; iMever et al 



2-4% for stars older than 1 Gyr (jBeichman et al 
20081 ). For A stars, a sim- 



ilar trend is seen, although the excess rate for 
young A stars (< 190 Myr) is even higher at 33% 



(Su et al 



2006). The age estimates for our targets 
stars (0 Leo: 50-380 Myr; £ Lep: 180-49 Myr; 
Lachaume et al. ( 19991 ); Chen et al. ( 20061 )) place 



them in or near the age brackets for the higher- 
percentage hot dust population. Some theoretical 
models of planet formation predict higher planet 
formation r ates for A stars, as compa red to solar 
mass stars ([Kennedy fc Kenvonlliooih . Thus, de- 
tection of hot dust in our targets systems is not 
necessarily unexpected. 

Other debris disks have near-infrared excesses 
detected through interferometry. Including this 
work, nine A and early F stars have been ob- 
served and, not including our tentative detection 
of £ Lep, three (Vega, £ Aq l and Leo ) were 
detected (ICiardi et all l200ll; I Absil et all l2006t 
di Folco et al.ll2007t I Absil et al.ll2008l) at the 1-2% 
excess level . Base d on limited radial velocity data, 
Absil et al. ( 20081 ) suggest that their detected ex- 
cess toward £ Aql may arise from a close M star 
companion and based on a reanalysis of the Spitzer 
data, they also find that there is no mid-infrared 
excess from this source. Only one lower mass 



star, r Ceti, has an observed near-infrared excess 
(|di Folco et al.ll2007t l. although the observations of 
FGK stars are currently limited by the sensitivity 
of the instruments. The near- infrared flux levels 
of the detected sources are all at the few percent 
level, although this may represent the brightest 
examples of a population of disks with hot grains, 
as the limits on the non-detected sources are not 



substantially lower, e.g. Idi Folco et al.l (120071 ) set 
an upper limit of 0.6% for the near-infrared emis- 
sion from e Eri. 

The dus t dist ribution inferred for Vega by 
Absil et al. (2006) is somewhat similar to Leo 



in that they modeled the inner dust with an in- 
ner radius of 0.17 to 0.3 AU with a very steep 
(n(r) oc r~ 4 ) radial power law, which is similar to 
a ring. The inner dust mass inferred was 8 x 10 -8 
with a dust luminosity of 5 x 10 -4 L Q . 
Although the number of debris disk stars sur- 
veyed is still relatively small, there are three detec- 
tions (one marginally significant) of near-infrared 
excess for which a population of hot dust is the 
most likely explanation, and here we consider if 
they represent a stable dust population produced 
by collisions of larger bodies orbiting close to the 
star or are the result of a transient event. The pro- 
duction of small, hot dust grains from the break- 
up of a comet or asteroid has been invoked in 
other cases, such as HD 69830 (IBeichman et al 



2005). Using a density of 3 gm cm -3 , a mass of 
5 x 10 -9 M ffi can be generated from the break-up 
of a single body with a 10 km radius. 

To evaluate the likelihood of so many tran- 
sient events in A stars, we used the model of 



Wvatt et al.l (|2007l ) for the steady-state evolution 



of collision-dominated debris disks, in which they 
derive a maximum possible dust luminosity as a 
function of age. In examining the properties of 
46 known A-star debris disks, they found only 4 
stars with a dust luminosity significantly higher 
than this maximum. A dust luminosity well above 
the steady-state maximum suggests either a tran- 
sient origin for the dust or unusual properties 
for the planetesimal belt. The mid-infrared ex- 
cess of £ Lep is above this threshold, while the 
mid-infrared excess of Leo is below. We use 
the same calculation to evaluate the inner dust. 
For simplicity, we take the maximum stea dy-state 
flux (/ max ) derived bv I Wvatt etaP (|2007l ) for the 
outer planetesimal ring and scale that value to 
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Model parameter (3 Leo C, Lep 



Inner ring 






r in (AU) 


R su b - 0.2 


R su b - 0.2 


Ar (AU) 


< 0.5 


< 0.15 


t(2/ito) 


2-20 xl0~ 3 


1-5 xl0~ 3 


a 


no constraint 


<-l 


X ((um) 


<2 


<2 


Outer ring 






n„ (AU) 


7- 15 


0.5 - 1.2 


Ar (AU) 


3 - 10 


11 - 15 


r(2/im) 


1 - 4 x 10~ 3 


2-3 xl0~ 4 


a 


no constraint 


-0.2 - 0.1 


A (A*m) 


35 - 70 


> 30 


5 


2 


2 



Table 8: The model parameters values from fitting the SED, IRS data and K-band excess. The ranges given 
cover a 67% probability range. 



the inner ring radius using their derived relation- 
ship / max oc r 7 / 3 . This ignores any difference in 
grain properties between the disks but is accept- 
able for an order of magnitude calculation. For 
both (3 Leo and £ Lep the near-infrared dust lu- 
minosity (/) compared to the maximum allowed is 
f //ma x ~ 10 6 , obvi o usly a bove the threshold of 10 
set bv lWvatt et al.l (|2007l) for anomalous systems. 
We also calculated this quantity for Vega, which 
for the inner dust also has a value of ///max ~ 10 6 . 
In this model of collisionally-dominated disks, all 
three of these objects are orders of magnitude 
higher than the expected steady-state flux, sug- 
gesting a transient event as the most likely origin. 
The finding from our simple dust model of small 
dust within the inner ring, but not the outer, mid- 
infrared producing ring, may also favor a transient 
origin for the near-infrared producing dust given 
the issues of dust lifetime to radiation pressure 
and sublimation. A recent planetcsimal collision 
or comet passing would drastically change the dust 
radii distribution and dynamics. 

Although more objects should be sampled to 
come to a stronger conclusion, it is suggestive 
that the near-infrared excess in these objects arises 
from a recent collision or cometary passing event. 
The near- infrared excess is an ideal probe of hot, 
small grains in these systems, as dust much closer 
to the star will sublimate. An observational test 
of the hypothesis that the near-infrared flux arises 
from emission from grains near the sublimation ra- 



dius is to make observations at other wavelengths, 
particularly H (1.6 /jm) and L (3.5/im) bands, to 
probe the wavelength dependence of the excess. If 
the flux is dominated by emission, the peak will be 
near K and L, while if it is scattering, the excess 
at H will be much higher and the excess at L much 
lower. 

6. Conclusions 

We have presented near-infrared interferometry 
observations of two A stars, (3 Leo and ( Lep, 
which were known to have mid-infrared excess 
emission from a debris disk. A near-infrared ex- 
cess of 1-2% was detected, although the detection 
for £ Lep should be confirmed. The interferometer 
observations do not spatially resolve the emission 
distribution, but place a maximum on the radial 
extent through the field of view, and in conjunc- 
tion with the spectral energy distribution, the spa- 
tial distribution of dust can be constrained. Both 
objects can be modeled as having a thin ring of 
dust grains at or near the sublimation radius in ad- 
dition to the previously known mid-infrared emit- 
ting belt. 

Although the models presented here are not a 
unique fit to the data, particularly with respect 
to the grain population, we can place strong con- 
straints on the dust composition and morphol- 
ogy. Both objects require small, non-silicate grains 
to be consistent with the near-infrared and mid- 
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Fig. 2. — The SED and flux excess for (3 Leo (top 2 panels) and £ Lep (bottom 2 panels). For each object, the 
top plot shows the Kurucz-Lejeune model used for the stellar photosphere as a solid line, with photometry 
from SIMBAD shown as points and the Spitzer IRS data as a thick line with errors. In some cases the error 
bars are smaller than the points. Our disk model is shown as a dashed line. The bottom plot for each object 
shows the disk model (solid line) with the CHARA, IRS, and MIPS excesses. For (3 Leo, the model shown 
has an inner ring with r- m = 0.13 AU, Ar = 0.3 AU, Tinner = 3.8 x 10~ 3 , thinner = — 1.5, A inner = 2/im and 
outer ring with r- ln = 13 AU, Ar = 6.2 AU,T ou tcr = 3.8 x 10~ 4 , ctWer = -1.5, A G uter= 50 /xm. For £ Lep, 
the model shown has an inner ring with r ln = 0.16 AU, Ar = 0.05 AU, T m — 4.5 x 10~ 3 , a; n =-2, A D ; n =2 /jm 
and outer ring with r in = 0.8 AU, Ar = 13 AU, r oute r = 2.6 x 10~ 4 , a outC r = -0.1, A D outer= 35 ^m. Both fits 
use a grain emissivity of e = e (A/A )~ 2 . 



infrared excesses. The minimum grain size re- 



quired (<~ 0.1 /im) is an order of magnitude smaller 



18 



than the nominal radiation pressure blowout ra- 
dius for spherical grains and requires a high pro- 
duction rate of small grains if some depletion 
does occur due to radiation pressure and subli- 
mation. For (j Leo, the near-infrared excess can 
not arise from dust generated by the planetesimal 
belt which produces the mid-infrared excess. For 
C Lep, the most likely model using simple geomet- 
ric distributions is also two separate rings, but it 
is possible that the larger bodies producing the in- 
ner dust may form a continuous extent with the 
mid-infrared planetesimal belt, although not with 
a simple, single power-law radial distribution. The 
luminosity of the inner dust is exceptionally high 
in comparison to steady-state evolutionary models 
of collisionally-dominated debris disks, suggesting 
origin in a transient event, such as the break-up 
of a comet or asteroid near the star. 

These observations are limited by the sensi- 
tivity of the current instruments, but with im- 
provements in near-infrared intcrferometry and 
other techniques, such as nulling intcrferometry 
and adaptive optics with coronagraphy, it should 
be possible to survey many more stars to deter- 
mine the population with hot dust. 
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